Justification of Power-Law Canonical Distributions 
Based on Generalized Central Limit Theorem 
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^ A self-consistent thermodynamic framework is presented for power-law canonical 

distributions based on the generalized central limit theorem by extending the discussion given by 
Khinchin for deriving Gibbsian canonical ensemble theory. The thermodynamic Legendre 
transform structure is invoked in establishing its connection to nonextensive statistical 
mechanics. 
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Systems exhibiting power-law behavior in their distributions are often encountered in nature. 

Until very recently, understanding of such properties based on the principles of Gibbsian 

canonical ensemble theory has not been satisfactory due to the exponential form of its 

distributions. Gibbsian canonical ensemble theory has been variously elaborated by several 

approaches such as the method of steepest descents, counting algorithm, maximum entropy 

principle, and the ordinary central limit theorem. (See Refs. [1,2] for a comprehensive account.) 

Among all these approaches, the one based on the central limit theorem seems to be the most 

satisfactory with the least number of underlying assumptions [3]. In this Letter, we present the 

derivation of the non-Gibbsian power-law canonical distributions based on the generalized 

central limit theorem of Levy and Gnedenko [4]. In the spirit of Khinchin [3], we thus establish 

the clear role played by probability theory in statistical mechanics of physical systems obeying 

power-law distributions. 

Suppose that a given system is composed of a large number of subsystems with energies 
{£, = EJB N > 0} , where B N is a positive Af-dependent factor which will be determined 

subsequently. The total energy is assumed to be the sum E = £ } + £ 2 +L +£ N . Here, it is useful to 
recall the ordinary central limit theorem. If each of E ( 's obeys a common distribution /(£,) with 
ordinary finite second moment (-E 1 , 2 ^ = | dE i E t 2 /(•£',■), then its iV-fold convolution, 

B N (f*f*L *f)(B N E), approaches the Gaussian distribution in the limit of large N, where 
{f*s)( x )=\ dx' f(x - x') g(x'). This property was exploited by Khinchin to establish 

Gibbsian canonical ensemble theory within the framework of probability theory. However, there 
are distributions for which all of the ordinary moments, (•E , i ") 1 = | q dE i E" /(is,-) (n = 1, 2, 3, L ), 

are divergent. In spite of this property, there still exists the limiting distribution, to which the N- 

fold convolution of the distribution of this type may converge [4]. Mathematically, this statement 

is called the generalized central limit theorem and the associated limiting distribution is referred 
to as the Levy-stable distribution, denoted here by F (E^). The latter is the case which leads to 

the power-law distributions of concern in the present work. The characteristic function of the 
Levy distribution in the half space is known to be of the form [4] : 
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Xr (t) = ]dE t e tE « F r (E) 

(1) 



= exp<^ -aid 7 exp 



'sgn(t) — 



where a is a positive constant, 7 the Levy index lying between (0, l), 9 a constant satisfying 
I 6 1 < 7, and sgn(?) = t /| t | the sign function of t. This is an invariant distribution in the sense 
that if each of the 2s ( .'s obeys i 7 r (£' / .), then their Af-fold convolution B N {F y *F y *~L *F y }(B N E) 
has the same characteristic function as that of the original F y (£ , / ), i.e., X^if) = [/Cy( ? /^v)] = 
X y {t), if B N is chosen to be 



B N = N l,r . (2) 



F y (E) is found to exhibit the following power-law behavior for large values of Er. 



F r (E t ) ~ E+y. (3) 



Now, we consider a power-law distribution 



5-1 



where is a positive constant depending on a certain parameter The range of s is 

determined in such a way that /(is,-) is normalizable and all of the ordinary positive integer 

moments are divergent. The logarithm of the characteristic function of the iV-fold convolution of 
this distribution with the same choice of B N as in eq. (2) is calculated to be 
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lnx <N) (t) = -a\t\" 'exp 



/sgn(f) 



On 



+{s-\)^\t\ s - l J dfy(exp[/ysgn(0]-l) 



y + 



B 



■n J 



■ (6) 



Here, a and are respectively given by 



a = (s -I) % s - l {P)\L(s -l)\>0, 



(7) 



= 1-5, 



(8) 



where [4] (p. 169) 



L («) = |^( e " y - 1 )<° (o<«<i). 

^ 



(9) 



The last term in eq. (6) is the correction term, which vanishes in the limit N — > °°. In deriving 
the above result, we have compared eq. (6) with eq. (1) and have found the Levy index to be 



7 = 5-1. 



(10) 



To completely characterize the limiting distribution F y (E), it is necessary to relate the 

constant a to a statistical physical quantity. Note that, in contrast to Khinchin's discussion based 
on the ordinary central limit theorem, here we have no finite ordinary moments of the 
distribution. This calls for a different approach to the determination of the constant a appearing 
in eq. (6). Our strategy is, therefore, to fix the constant £ by some means using the distribution 
of a single subsystem. 
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To this end, we observe that although the ordinary first moment, (£',■) , is divergent the 
generalized first moment, (E^ = | dE i E l P q {E), defined in terms of the escort distribution [5] 

associated with /(£,), that is, 

f(E) 

J ^ 



is designed to be finite if q is chosen to be 



4 = 1 + -. (12) 

s 



This choice is dictated by the minimal requirement from the finiteness of {E^ . Let us calculate 
the denominator in eq. (1 1) as well as (E t ) : 

c q ^]dE i r(E i )= { -^^"(p), (13) 



Therefore, we clearly see how the constant a in eq. (7) is expressed in terms of u , which is a 

physical quantity, namely the generalized internal energy. From eqs. (13) and (14), it also 
follows that 



', = —<"*■ (15) 



In order to develop the connection with thermodynamics, we regard /? and u q as a conjugate pair 

in the thermodynamic Legendre transform structure [6]. To do this, we define an entropy 
functional, S q [f], which satisfies 
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¥--P. (16) 



dT 

W su " <17) 



where 



r q = s q -p Uq . (18) 

This requirement determines the relation between £ and f5 appearing in eq. (4), once the form of 
the entropy S [f] is specified. At this juncture, we examine the following form 



SJ/^T^-K^- 1 )' ( 19 ) 



where we have introduced a scale o to make S q [f] dimensionless. Then, using eq. (16), we 
obtain 

P = — ^h- (20) 

From eqs. (18) and (20), we verify that equation (17) indeed holds, as it should do. 

We complete our discussion by pointing out that actually the entropy in eq. (19) and the 
associated distribution of a single subsystem in eq. (4) are respectively the Tsallis nonextensive 
entropy [7] and its optimal distribution under the constraint on the generalized first moment with 
respect to the escort distribution in eq. (11) [8]. We mention that similar studies on convergence 
properties in both the full and half spaces have recently been reported in the context of 
nonextensive statistical mechanics in Refs. [9,10]. 
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In conclusion, we have constructed a self-consistent thermodynamic framework for power- 
law canonical distributions by making use of the generalized central limit theorem in a similar 
manner to that given by Khinchin for deriving Gibbsian canonical ensemble theory. In 
establishing a connection between generalized canonical ensembles with power-law distributions 
and nonextensive statistical mechanics, we have utilized the thermodynamic Legendre transform 
structure, which is an underlying principle common in both extensive and nonextensive 
statistical mechanics. 
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